Coherent Logic Gate for Light Pulses based on Storage in a Bose-Einstein Condensate 
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A classical logic gate connecting input and output light pulses is demonstrated. The gate operation 
is based on three steps: First, two incoming light pulses are stored in a Bose-Einstein condensate, 
second, atomic four-wave mixing generates a new matter wave, and third, the light pulses are 
retrieved. In the presence of the new matter wave, the retrieval generates a new optical wave. The 
latter will only be generated if both input light pulses are applied, thus realizing an AND gate. 
Finally, we show that the gate operation is phase coherent, an essential prerequisite for a quantum 
logic gate. 



Single photons are well suited for quantum commu- 
nication over long distances. To perform quantum in- 
formation processing with single photons, however, one 
must find a physical process in which a single photon 
drastically alters some property of another single photon. 
This is a major challenge because in traditional nonlin- 
ear optical media the nonlinearities are much too weak to 
generate an appreciable effect on the single-photon level. 
Several techniques for addressing this problem have been 
proposed and are being pursued experimentally, namely 
the use of atoms in optical resonators the use of 

additional light to drive Raman transitions in atoms 0- 
5], and the use of the dipole-dipole interaction between 
Rydberg atoms 

Here we present a first experiment that explores the 
avenue of generating a logic gate for classical light pulses 
by temporarily converting the light pulses into atomic ex- 
citations in a Bose-Einstein condensate (BEC) and using 
s-wave collisions between pairs of ground-state atoms. 
These collisions are responsible for the appearance of 
the nonlinear term in the Gross-Pitaevskii (GP) equa- 
tion. In the context of quantum information processing, 
they have been used to generate massive entanglement 
between many atoms @] but not to generate a logic gate 
for light pulses. In addition, we demonstrate that the 
gate operation is phase coherent, an essential prerequi- 
site for a quantum logic gate. 

We use a geometry in which the nonlinearity of the GP 
equation creates a new atomic momentum component 
by four-wave mixing (FWM) of matter waves [10Hl3j in- 
volving two spin states [3 GJ]. Upon mapping the new 
atomic momentum component back onto light, it creates 
population in a new optical momentum component. This 
light emission process is accompanied by Raman amplifi- 
cation of matter waves (AMW) [II Hz)]. The light emit- 
ted during Raman AMW and the phase coherence of this 
light have never been studied experimentally; despite re- 
lated work in atomic FWM []J], Rayleigh AMW |ml2Cj. 
and superradiant light scattering [2lj, [22j . 

A scheme of our experiment is shown in Fig. [T] The 
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87 Rb hyperfine states |1) = \F = 1) and |2) = \F = 2) 
of the 55i/2 ground state together with the |e) = \5P\i2) 
excited state, each with rap = — 1, form a A scheme in 
which Raman transitions are driven. More precisely, the 
Raman light fields are tuned to the two-photon resonance 
with the single-photon detuning chosen exactly midway 
between the |1') = \F' = 1) and \2') = \F' = 2) compo- 
nents of the excited state |e). 

Initially, a BEC with N ~ 1.5 x 10 6 atoms is prepared 
in an optical dipole trap with measured trap frequencies 
of (w x ,u y ,u z ) — 2tt x (70,20,20) Hz, in internal state 
|1), and at zero momentum, which we denote as |1, 0). A 
magnetic hold field of ^ 1 G is applied along the z-axis 
(orthogonal to the plane shown in Fig. [IJ to preserve the 
atomic spin orientation. All light fields applied in the 
experiment drive tt transitions. 

After preparation, the BEC is illuminated by a Raman 
pulse, consisting of signal light with wave vector fc s , prop- 
agating rightward in Fig. [TJ and control light with wave 
vector fc c i, propagating upward. Signal light is absorbed 
and coherently stored in the atomic state |2,fc) with in- 
ternal state 1 2) and wave vector k = k s — fe cl . The pulse 
area of the Raman pulse is chosen such that ~ 1/3 of the 
atomic population is transferred to state |2, k). 

Immediately thereafter, a second Raman pulse is ap- 
plied with the direction of the signal light as before, 
but now with the control light propagating leftward with 
wave vector k C 2- During this pulse, signal photons are ab- 
sorbed and stored in state |2, k + q) with q = k s k C 2 — k 
23]. This Raman pulse has a duration of ~ 100 fis. The 
pulse area of ~ 7r/2 yields equal populations of states 
|1,0), |2,fc), and |2,fc + g). 

In principle, the second Raman pulse could simulta- 
neously drive a second process in which population is 
transferred from state |2, k) to state |1,— q). In prac- 
tice, however, the nonzero initial momentum creates a 
Doppler shift for the resonance frequency of this process. 
By fine tuning the two-photon detuning we resonantly 
drive the |1,0) — > |2,fe + q) processes and drastically 
suppress the |2, k) —> \l,—q) process. The pulse is long 
enough to make interaction-time broadening small com- 
pared to the energy splitting between the two resonances 
frequencies. 

During the following dark time, with duration £fwm, 
atomic FWM with two internal states populates the state 
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FIG. 1: (Color online) Scheme of the experimental pro- 
cedure. Two Raman pulses are applied, preparing three 
BECs with momenta 0, fc, and k + q. During a subsequent 
dark time, atomic FWM creates a BEC with momentum q. 
Arrows show the propagation directions of the light beams. 
Circles represent atomic momentum components. The inter- 
nal state is color coded: white=|l), red=|2). A third light 
pulse with only control light applied retrieves the signal light. 
The retrieved light will have a component propagating down- 
ward only if signal light is applied during both Raman pulses, 
thus realizing an AND gate. A modified state preparation can 
additionally populate the momentum component — q so that 
atomic FWM also populates fc — q. This extended scheme, 
in which the dashed circles are populated, generates a time- 
dependent interference pattern in the retrieved light. 



1 1,9). The FWM can be understood intuitively as an 
atomic scattering process. An atom in state |1,0) col- 
lides with an atom in state |2, k + q). The existing BEC 
in state \2, k) creates bosonic enhancement for one atom 
to emerge in this state. Conservation of momentum and 
of the internal-state energy makes the other atom appear 
in state 1 1 , qr) . In addition, conservation of kinetic energy 
requires k ■ q ~ [HI [l2| . Note that the FWM process 
in our experiment is clearly distinguishable from spin ex- 
change, unlike the only previous experiment on atomic 
FWM with two internal states [HJ]. The Bose-enhanced 
creation of other atomic momentum components \l,nq) 
and 1 2, k + nq) with integer n would conserve momen- 
tum but not kinetic energy and is therefore negligible. 
The FWM occurs inside the optical dipole trap in order 
to avoid a slowdown of the FWM due to the reduced 
density in a mean-field driven expansion. 

To test whether atomic FWM occurs in our experi- 
ment, we study time-of-flight absorption images. The 
inset in Fig. [2] shows such an image. Four atomic momen- 
tum components are clearly distinguishable. The addi- 
tional tiny signal in the bottom left corner shows that the 
process populating state |1,— q) during the second Ra- 
man pulse is suppressed drastically but not completely. 

To confirm that the atom number N q in the momen- 
tum component \l,q) is actually generated by FWM, 
we study its temporal growth. To avoid FWM during 
the time of flight, we abort the FWM by applying an 
8 /t/s pulse of depletion light [HI, 65 MHz blue detuned 
from the |2) <H- |1') transition, immediately before release 




FIG. 2: (Color online) Atomic four-wave mixing. The 

atom number N q in the atomic momentum component q is 
extracted from time-of-flight absorption images. N q depends 
on the four-wave mixing time. At short times, this depen- 
dence is quadratic. Inset: An image oriented as Fig. [T] taken 
for £fwm = 1.8 ms without applying depletion light. 



from the dipole trap. The atom number N q extracted 
from time-of-flight images taken after applying depletion 
light is shown in Fig. [5] For short times, N q displays a 
quadratic growth with time, as expected for FWM. For 
longer times, N q saturates because the different momen- 
tum components no longer overlap spatially. The line 
shows a fit to the data, where the observed timescales 
for the initial growth and for the saturation agree fairly 
well with theory [23[ . 

After confirming that N q is actually generated by 
FWM, we now map the atomic states back onto the light 
field. To this end, the BEC in the trap is illuminated 
by a third light pulse, the retrieval pulse, during which 
a control beam propagating upward is applied, whereas 
no signal light is applied. We use a detuning of 300 MHz 
red from the |2) f-> |2') transition because we find ex- 
perimentally that this maximizes the retrieved photon 
number propagating downward. Using control light with 
an intensity of roughly 100 mW/cm 2 , we find that the 
retrieved light emerges in less than 5 /is. Each atom in 
internal state |2) is transferred back into internal state 
|1) in a Raman process under emission of a signal pho- 
ton. These Raman processes are bosonically stimulated 
by the two BECs in states |1,0) and |l,<jf). The stimu- 
lated growth of atomic population in these two BECs is 
called Raman AMW. Along with the population growth 
in states |1,0) and 1 1 , qr) , signal light with two momen- 
tum components is emitted, one propagating downward, 
the other rightward. 

We concentrate on the downward propagating com- 
ponent with wave vector k s — q. It is created when 
atoms are transferred from state |2, k) to |1, q). This will 
only be possible if atomic FWM occurs because otherwise 
N q = 0. FWM, in turn, will occur only if the signal light 
is applied during both Raman pulses. We use an electron- 
multiplying charge-coupled device (EMCCD) camera to 
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FIG. 3: (Color online) Logic gate. The retrieved photon 
number propagating downward is shown for four different ex- 
perimental settings, in which the signal beams during Raman 
pulses 1 and 2 are turned on or off independently. Downward 
propagating light will only be retrieved if both signal beams 
are on. This demonstrates an AND gate for classical light 
pulses based on storage in a BEC and FWM of matter waves. 



measure the photon number propagating downward. An 
iris diaphragm is placed in an intermediate imaging plane 
to suppress stray light from the control beam. Results 
for £fwm = 0.4 ms are shown in Fig. [3] They clearly 
demonstrate an AND gate for the two classical signal 
light pulses. 

The rightward propagating component of the retrieved 
light is of less interest here. Its existence does not rely 
on atomic FWM. If the second pulse were omitted com- 
pletely, there would still be the retrieval of rightward 
propagating light, well known from experiments on elec- 
tromagnetically induced transparency (EIT) [24J. The 
photon number retrieved in this beam, however, does de- 
pend on whether FWM occurred because there is compe- 
tition between the processes retrieving light propagating 
downward and rightward. 

When considering the perspectives for scaling this gate 
down to the single-photon level in order to obtain a quan- 
tum logic gate, it is crucial whether the gate operation is 
phase coherent. We will now show experimentally that 
this is the case for the gate demonstrated here. 

From a theoretical point of view [I?} , Raman AMW is 
analogous to a usual stimulated Raman process, except 
that the emission is bosonically stimulated not by appli- 
cation of a second laser beam but by the presence of a 
second BEC. Effectively the role of the second laser beam 
and the second BEC are exchanged. In a usual stimu- 
lated Raman process, the atoms are transferred into the 
initially empty state in a phase coherent way, with the 
relative phase of the two applied laser beams determining 
the phase of the transferred atomic amplitude. Similarly, 
we expect that Raman AMW generates light in a phase 
coherent way, with the relative phase of the two BECs 
determining the phase of the emitted light. 




FIG. 4: (Color online) Phase coherence. The retrieved 
photon number propagating downward exhibits a sinusoidal 
oscillation as a function of the FWM time. The best-fit value 
for the oscillation frequency is 15.4 kHz. The oscillation oc- 
curs because two pathways contribute to this signal in the 
extended scheme with six momentum components. 



To test experimentally whether the emitted light is 
phase coherent, we shorten the second Raman pulse 
to ~ 35 /zs at correspondingly higher light intensities. 
Interaction-time broadening now makes the Doppler shift 
for the transfer of population into state |1, —q) irrelevant. 
We choose pulse areas of ~ tt/2 for both Raman pulses to 
create four equally populated atomic momentum compo- 
nents 1 1 , 0) , \l,—q), \2,k), and \2,k + q). Subsequently, 
atomic FWM populates the states \l,q) and \2,k — q), 
see Fig. [TJ We refer to this as the extended scheme be- 
cause it features six atomic momentum components, in- 
stead of four. The retrieval pulse, applied as before, again 
will generate downward propagating signal light only if 
atomic FWM occurs. But now, two pathways contribute 
to this signal. The light can be generated by transfer of 
an atom either from state |2, k) to |1, q) or from |2, k — q) 
to|l,0). 

The retrieved photon number propagating downward is 
shown in Fig. 0] as a function of £fwm- It clearly shows a 
sinusoidal oscillation. To understand the physical origin 
of this oscillation, we note that during £fwm, the atomic 
components that contribute to the first pathway, |2, fe) — > 
\l,q), differ in kinetic energy by AE\ = h 2 (k 2 — q 2 )/2m, 
whereas the components of the second pathway, |2, k — 
q) -> |1,0), differ by AE 2 = Ti 2 {k~ q) 2 /2m. Here m 
denotes the atomic mass. During £fwm, these pairs of 
BECs thus accumulate different relative phases. Upon 
retrieval, these relative phases of the pairs of BECs are 
mapped onto the phases of the retrieved signal light, as 
discussed above. Eventually, the amplitudes associated 
with the two pathways create an interference signal on 
the detector. An independent measurement of our beam 
geometry yields q 2 = 2.08k 2 and q ■ k = -0.037k 2 . With 
2ir/k s = 794.979 nm we expect an angular frequency u> = 
(AE 2 -AEi)/H = 2?rx 15.4 kHz for the oscillation. If the 
beams pointed exactly along the coordinate axes, then we 
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would expect oj — 4E lec /h, where E Iec — h 2 k^/2m is the 
recoil energy of the signal light. 

Two effects contribute to the observed envelope of the 
oscillation 23]. First, Raman AMW requires spatial 
overlap of the atomic momentum components. For long 
times, such overlap is lost and no directed retrieval is ob- 
tained. The timescale seen for this effect in Fig. U] is sim- 
ilar to the timescale of saturation in Fig. [31 as expected. 
Second, retrieval of downward propagating light requires 
FWM to occur, so that an increase at short times, as 
seen in Fig. [31 is also seen in the envelope in Fig. [H 

For simplicity, we fit a sinusoid with a Gaussian en- 
velope and a constant visibility V to the data in Fig. [4] 
This yields best-fit values of uj/2tt = 15.4 ± 0.1 kHz and 
V = 0.35 ± 0.02. The fact that V is not perfect is the- 
oretically expected 23]. The agreement with the above 
expectation for w is excellent. This oscillation is not seen 
in time-of-flight absorption images, either in Fig. [2] or in 
similar data that we took for the extended scheme (not 
shown here) . This demonstrates that the retrieval of the 
light is necessary to make the oscillation appear. 

The observed oscillation proves that the complete gate 
operation is phase coherent. This includes all physical 
processes involved in the gate, namely Raman pulses 1 
and 2, atomic FWM, and light emission during Raman 
AMW. The observed phase coherence is a crucial ingredi- 
ent for a possible extension to a photon-photon quantum 



logic gate. 

It should be noted that an extension of our scheme to 
the single-photon level is challenging. One possible prob- 
lem is that smaller photon numbers slow down the atomic 
FWM so that the spatial overlap of the matter waves 
might end before achieving enough population transfer. 
Another possible problem is that the bosonic stimulation 
caused by a single atom in driving either atomic FWM or 
Raman AMW is weak so that competing processes might 
cause problems, e.g., isotropic s-wave scattering or spon- 
taneous photon emission. Working out strategies to solve 
these problems is beyond the scope of the present paper. 

On a more general level, our experiment demonstrates 
that storage and retrieval of light combined with the 
nonlincarity of the GP equation can be used to cre- 
ate a phase-coherent gate for two classical light pulses. 
If the geometry were altered to make all laser beams 
copropagating, then the nonlinearity of the GP equa- 
tion would generate a conditional phase shift instead of 
atomic FWM. A recent theoretical analysis of the single- 
photon version of that scheme resulted in a detailed pro- 
posal for a photon-photon quantum-logic gate (25j. 
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APPENDIX 

I. DYNAMICS OF FOUR- WAVE MIXING 

A. Initial Growth of N q (t) 

The growth of N q (t) at short times can be estimated as 
follows. Let No, N q , Nk, and N k+q denote the atom num- 
bers with the various momenta, a\2 the interspecies scat- 
tering length, and g\ 2 = A / KTb 2 ai 2 /m. We write the single- 
particle wave-function of a Thomas-Fermi parabola as 



1/2 



(1) 



where this is real and u(x) = otherwise. The normal- 
ization condition J d 3 x\u(x)\ 2 = 1 yields 



V — —R x RyR z 



(2) 



V has the dimension of a volume and expresses the ratio 
of particle number over peak density. From the measured 
trap frequencies and the measured initial total particle 
number, the Thomas-Fermi radii in our experiment are 
estimated to be (R x , R y , R z ) — (8, 27, 27) /im. 

Based on coupled Gross-Pitaevskii equations [SlJ , one 
can generalize Eq. (17) of Ref. [H| to the case of FWM 
with two internal states. This yields the following result 
for the condensate wave function of the component 1 1 , q) 
at short times 

^ q (x,t) = ^u(x)\u(x)\ 2 ^N (0)N k (0)N k+q (0). (3) 

This is normalized such that N q (t) = J d 3 x\ip q (x, t)\ 2 . 
Using 



d 3 x\u(x)\ ( 



21V 2 



(4) 



we obtain 

where we introduced the timescale for the initial growth 
hV I 



21 



g 12 V 8N (Q)N k (0)N k+g (oy 



(6) 



For [iV (0)iVfe(0)7V fe+? (0)] 1 / 3 = 2.8 x 10 5 and a l2 = 98.4 
Bohr radii [Sj|, Eq. © predicts r = 2.1 (is. 



B. Fitting Function for Fig. [2] 

At longer times, FWM can make the atomic popula- 
tion oscillate back and forth between different momentum 



components, at least in principle S2]. In our experiment, 
however, the FWM dynamics come to an end because 
the different momentum components spatially separate 
as time progresses. The FWM dynamics is terminated at 
a time when the fraction of the atomic population that 
has been transferred is still small and no oscillations have 
occurred yet. 

Using two coarse approximations, we now derive an 
analytic expression for N q (t) for all times. To represent 
the decreasing spatial overlap of the clouds at long times, 
we take the different center-of-mass momenta of the mo- 
mentum components into account. However, we neglect 
a possible time dependence of shape, radii, and particle 
numbers of the three initially populated Thomas-Fermi 
parabolae. Hence, the time dependence of the conden- 
sate wave functions for a £ {0, k, k + q} is approximated 
as 



hk n 



(7) 



with k a = 0, k, k + q for a = 0, k, k + q, respectively. 
In our experiment, the FWM transfers only a small frac- 
tion of the population, all relevant scattering lengths are 
almost identical, and the timescale for the clouds to sep- 
arate spatially is much shorter than the trapping period. 
Hence these approximations seem reasonable. 



Furthermore, we simplify Eq. (15) of Ref. [S2| by ne- 
glecting mean-field energies, external potential, and dis- 
persion of the newly generated wave packet. This yields 



d t + — ■ v 

m 



i>q = -^-Mk^k+q- (8) 



This is our first coarse approximation. It drastically sim- 
plifies the problem. For g i2 — 0, the solution would 
be a classical drift of the initial wave packet tp q (x,t) = 
ip q (x — hqt/m,0), much like Eq. (JTJ). For g 12 ^ 0, each 
drift trajectory x(t) = x(0) + hqt/m accumulates an am- 
plitude due to FWM but there is no crosstalk between 
different trajectories because we neglected the dispersion 
of the wave packet. The initial condition is ip q (x,Q) = 0. 

Considering the trajectory that originates at the cloud 
center £e(0) = 0, we obtain 



hq 



f.t 



912 



dt 3 ipo(x 3 , t 3 )ipl(x 3 , t 3 )ip k+q (x 3 , t 3 ) (9) 



with x 3 = hqt 3 /m. We approximate the laser beam 
geometry as rectangular with 



1 



(10) 
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This yields 



II. FRINGE VISIBILITY AND ENVELOPE 



ip q — t,t 




*2(t) 



du 1 




(11) 



where we introduced two versions of the timescale which 
describes saturation 



to 



mRy 
2hk, 



2t 



(12) 



to and t\ are connected by the aspect ratio e = R y /R x . 
The measured trap frequencies yield e = Lo x /Lu y = 
70 Hz/20 Hz = 3.5. Eq. O predicts t = 2.3 ms. 

The condition that Eq. (fT]) will be used only if it is 
real, translates to the upper bound of the integral 



t 2 (t) = min{t,t ,ti}. 
The integral has the analytic solution 



1 TO / IT 




(13) 



(14) 



with 



f(x) = l( 2 ;v / r^[17 + e 2 -2.T 2 (l + e 2 )] 



+ (15 — e 2 ) arcsinx 



(15) 



Rather than performing a similar calculation for each tra- 
jectory and eventually carrying out a spatial integral, we 
assume that this trajectory is sufficiently representative 
for the whole wave packet. This yields our second coarse 
approximation 



N q {t) oc 



( —t,t 
1 TO 



(16) 



The factor of proportionality is obtained by matching 
the behavior at short times with Eq. Using f(x) = 
x + 0(x 3 ), this yields 



N q (t) 



f 



h(t) 



to 



(17) 



Finally, we add an empiric overall offset Ar o ff S ct to this 
expression for N q (t). Hence, our model has three free fit 
parameters: t, to, and N s sct . 

A fit of this model to the data in Fig. [5] yields the 
best-fit values r = 2.5 /is and to — 1.3 ms. We find 
good agreement with the prediction r = 2.1 /is. The 
coarse approximations used to derive the fitting function 
are most likely the reason why the agreement with the 
prediction to = 2.3 ms is not better. 

Both fit parameters have a simple graphical interpreta- 
tion in Fig. [2] First, the fit curve reaches N q (i) — N q (0) = 
10 4 approximately after the time t = r-y/10 4 = 0.25 ms. 
Second, the fit curve is time independent for t > t\ = 
2t /Vl + e 2 = 0.69 ms. 



A. Qualitative Discussion 

The fringe visibility in Fig. 0]is expected to be incom- 
plete. To understand this on a qualitative level, consider 
a hypothetical scenario in which control light during re- 
trieval would be applied for a very short time such that 
only a small fraction of the light would be retrieved and 
the atomic populations would be essentially unchanged. 
In this case, an appropriate choice of the initial atomic 
populations could easily balance the intensities emitted 
along the two pathways generating downward propagat- 
ing light. In this hypothetical scenario, one would ide- 
ally expect unit fringe visibility. If such a measurement 
yielded a fringe visibility much below unity, then one 
might be able to draw conclusions about the maximum fi- 
delity of a possible extension of the scheme to a quantum- 
logic gate. 

In our experiment, however, we apply control light 
during retrieval until no more signal light emerges be- 
cause this strongly improves the signal-to-noise ratio of 
the data and it still suffices to demonstrate phase co- 
herence in general. As a result, the atomic populations 
change drastically during retrieval. Specifically, at the 
end of the retrieval, there are essentially no atoms left 
in internal state |2). Neglecting spontaneous emission, 
this implies that all atoms originally in state |2, k + q) 
were transferred into state |1, q), which had a small pop- 
ulation before retrieval. In other words, the population 
of state \l,q) grows by a large factor during retrieval. 
The population growth in state |1,0), however, is given 
by a factor of ~ 2 at best, because this state already 
has a large population before retrieval and its growth 
comes dominantly from atoms in state |2, fc). Hence, the 
bosonic stimulation factors along the two pathways that 
contribute in Fig. U have very different time evolutions. 
Therefore, the intensities emitted along the two path- 
ways generating downward propagating light cannot be 
balanced throughout the whole retrieval process. This 
reduces the fringe visibility in the time-integrated sig- 
nal substantially. Hence, the fringe visibility measured 
in Fig. 0] does not allow for a trivial extrapolation for 
the fidelity of a possible extension of our scheme to a 
quantum-logic gate. 



B. Plane- Wave Model 



To show that we have a reasonable understanding of 
the value of the fringe visibility V measured in Fig. 21 
we now develop a simple model. We start by adapting 
the Hamiltonian in Eq. (7) of Ref. [s"ij to our geometry, 
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obtaining 

oo 

H = ^ h(Wb,J>lj>n + U c ,n4iCn) ~ ^ MjQcLj 
ri= — oo je{d,r} 
oo 

+ J2 (i%(4^+i^ + 44U+H.c). (18) 

n— — oo 

The operators b^ n and & n create bosonic atoms in states 
1 2, As + nq) and \l,nq), respectively. The first terms in 
H represent the kinetic energies of these atoms. The 
internal-state energies for atoms at rest are absorbed 
in an interaction picture. Equation (fTO)) yields Wb,„ = 
2(n 2 + l)u r and uj Cin — 2n 2 uj r , where uj r = E lec /h — 
hk 2 /2m. 

The operators d d and a\, create signal photons in plane- 
wave modes propagating downward and rightward, re- 
spectively. 5d and 5 r denote the two-photon detunings of 
the Raman transition. If we approximate all BECs and 
the control light applied during retrieval as plane waves, 
then the two emission modes will be fixed by momen- 
tum conservation. The terms in H containing Sd and S r 
express the photon energies in the interaction picture. 

The terms a^cjj +1 fe n and ajcjj6„ describe the emis- 
sion of signal photons propagating downward and right- 
ward, respectively. These processes are accompanied by 
the transfer of atoms from state |2,fc + nq) to states 
|1, (n+l)q) and |1, nq), respectively. The Hermitian con- 
jugate terms describe reabsorption of previously emitted 
signal photons. The coupling constant g is proportional 
to the Rabi frequency of the control laser applied dur- 
ing retrieval. Obviously, it suffices to consider g as real 
because the phase of g can be absorbed by resetting the 
phases of the operators ad and a r . 

We now consider the dynamics on a mean-field level. 
To this end, we calculate the Heisenberg equations of 
motion for the operators &j, b n , and c n . We take the 
expectation values of these equations, and obtain equa- 
tions of motion for mean fields Oj(i) = iflj) etc. We 
assume that the expectation values of products of oper- 
ators factorize, i.e. (cjj +1 & ra ) = (c„, i){b n ) etc. In addi- 
tion, we make a transition to another interaction picture 
with ctj = \/Naj exp(i5jt), b n — \/Nb n exp(— W(, >n t), 
and c„ = \f~Nc n exp(— ioj CiTl t), where N is the total atom 
number. We obtain 

d t a d = SivEX+Ae-^"*, (19a) 
d f a r = <7ivE„<&ne-^»*, (igb) 

d f b n = -g N (a d c n+1 e lAd nt + a r c n e iAr ' nt ) , (19c) 
d t c n = 5w (a^„_ 1 e- 4A ' i — lt + a ;fe„e- 2A ''-*)(19d) 

with g N = gVN, A d , n = 5 d + cJ&, n - w c ,n+i = Sd- Anui r , 
and A r ,„ = S r + uj]>. n — u) c ,n = $r + 2u r . 

In our experiment, g^ ^> tu r and all light is retrieved 
in a time very short compared to l/w r . Combined with 
the facts that \6jt\ <C 1 and that significant population 



in atomic modes exists only for small n, we conclude that 
|Ad >n i| <C 1 and |A r , n i| <C 1 for all relevant n. In the 
following, we therefore approximate all phase factors in 
Eq. (pi]) as unity. 

At this point, the model reproduces our qualitative ar- 
gument from Sec. Ill Al Consider a hypothetical scenario 
in which control light during retrieval would be applied 
for a very short time such that only a small fraction of 
the light would be retrieved and the atomic populations 
b n and c„ would be essentially unchanged. Specifically, 
we could ignore the build-up of population in atomic mo- 
mentum components with \n\ > 1. Here, the build-up of 
the amplitude ad would be given by 

d t a d = S-jv(cS&-i + c*6 )- (20) 

This equation expresses the quantum interference be- 
tween the two pathways. In this scenario, one would 
expect that an appropriate choice of the atomic ampli- 
tudes should ideally yield perfect fringe visibility. As ex- 
plained in Sec. Ill Al our experiment is not performed in 
this regime of a short retrieval pulse so that the tempo- 
ral evolution of the atomic populations b n and c n during 
retrieval needs to be taken into account. 



C. Irreversibility 

In our experiment, the BEC has a finite spatial exten- 
sion. As a result, emitted signal light can spatially leave 
the BEC, making a subsequent reabsorption impossible. 
This brings an effective irreversibility into the physical 
process which is not captured by the above plane-wave 
model. A full numerical model taking spatial wave pack- 
ets into account is beyond the scope of our present work. 
Instead, we add damping coefficients jd and j r to model 
this irreversibility. We replace Eq. (TT9")) by 

dtdd = -\ldO-d + ffwE„ c r*+l fo n: ( 21a ) 

d t a r = -\-i r a r + gNY< n c *n b n, (21b) 
d t b n = -gN {a-dCn+i + a r c n ) , (21c) 
dtc n = 9n {a*d b n-i + K b n) ■ (21d) 

The total numbers of photons Nd(t) and N r (t) irre- 
versibly emitted between times and t in each direction 
are easily calculated using 

d t Nd = N ld \a d \\ d t N r =N lr \a r \ 2 . (21e) 

The initial conditions are ad = a r — Nd = N r = 
and b n = c n — for |n| > 1. The times needed to 
travel through the BEC are proportional to the Thomas- 
Fermi radii in the corresponding direction, so that the 
ratio of the damping coefficients is determined by the 
inverse ratio of the corresponding Thomas-Fermi radii 
Ir/ld = Ry/Rx = e = 3.5. In general, the nonlinear 
differential equations (|2ip can produce fairly complicated 
non-harmonic oscillations. 
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As discussed in Sec. IIII A[ g = 2n x 0.6 MHz for our 
experiment. Combination with N ~ 10 6 yields g^ — 
2tt x 0.6 GHz. 

In the following, we use the simple estimate 7d = 
c/Ry ~ 1 x 10 13 s _1 with the vacuum speed of light c. 
Note that slow light (S5|, i.e. a substantial reduction of 
the group velocity of the signal light due to the presence 
of the control light, would only be obtained if the popu- 
lation of state 1 2) were small, which is not the case in our 
experiment. We obtain "fd/gN ~ 3x 10 3 . This is deeply in 
the regime jd/gN 3> 1, where the oscillations are over- 
damped. Adiabatic elimination of ad and a r from the 
equations (by formally setting dtad = dta r = 0) shows 
that the typical timescale for light to leave the BEC is 
Id I 'gff ~ 0-7 M s - This value agrees fairly well with our ex- 
perimentally observed timescale of ~ 1 /zs. In this regime, 
the value of "/d/gN only determines this timescale but has 
no other effect on the retrieval dynamics. Hence, even if 
the group velocity of the signal light were reduced, this 
would not affect the value of the visibility V predicted 
by the model, as long as one would stay in the regime 

id/gN > l. 

The initial conditions are determined by the pulse ar- 
eas $i and $2 of the two Raman pulses. The state before 
the first Raman pulse is described by c = 1. We assume 
that the first Raman pulse is so short that the kinetic en- 
ergies can be neglected during this pulse duration. The 
first Raman pulse then converts the state into 



Ci Si\(l 
Si Ci ) \ 



Ci 
Si 



(22) 



Similarly, the second Raman pulse converts this state into 



C-l 

b 



C2 —S2 

5*2 C2 

C2 — 02 
S2 C2 





(c x c 2 


»')= 


\ C1S2 




—S1S2 


in 


S1C2 



, (23a) 



(23b) 



cos % and Sk 



sin for 



Here we abbreviated CV- 

k G {1,2}. The subsequent FWM builds up population 
in b—i and c\. In addition, each state accumulates a 
phase proportional to its kinetic energy. This yields 



6-1 = -«/3cgC_i6o, 
& = e lx SiC 2 , 



bi 



C102, 



c_! = -e^Sift, 

Co = C1C2, 

ci = -ipb^bxCQ 



(24a) 
(24b) 
(24c) 



with x = — 2o; r iFWM and with a dimensionless parameter 
(3 that describes how far the FWM has evolved. For short 
i FWM , Eq. © yields /3 = y/8/2lg 12 t FWM N /hV. These 
initial conditions neglect the depletion of the other states 
caused by FWM because in our experiment |6-i| 2 <C 1 
and |ci| 2 < 1. 

We numerically solve the differential equations (f2"Tj) 
with the initial conditions (|24]l. To this end, we intro- 
duce a momentum cutoff by setting b n (t) — c n (t) = for 
|n| > 6. In addition to J r /^fd = 3-5 and jd/gN = 3 x 10 3 , 



we choose #1 = -d 2 — tt/2 and (3 = y/2. This yields 
|ci(0)| 2 = 1/32 ~ 0.03 which is realistic. The numerical 
calculation shows that some population is transferred to 
states with \n\ — 2, but population in states with |n| > 2 
is negligible. The calculation of Nd at long times yields 
an interference pattern as a function of x- The pattern 
is to a good approximation sinusoidal. A sinusoidal fit 
yields a best-fit value of V = 0.61 for the visibility. 

V is fairly sensitive to varying the parameters, e.g. 
V = 0.36 is obtained for (tfi,tf 2 ,/3) = (0.35?r, 0.65tt, y/2) 
as well as for (tfi,i? 2 ,/3) = (7r/2,7r/2, 0.3^2). The ap- 
proximations used in the model as well as the fact that 
the experimental control over the parameters is not per- 
fect let the experimentally observed value V = 0.35 from 
Fig. [4] seem reasonable. In particular, the finite visibility 
in this particular measurement does not indicate exper- 
imental problems that would force us to conclude that 
the fidelity of a possible extension of our scheme to a 
quantum logic gate would be limited. 



D. Envelope of the Interference Pattern 

Now, we develop a simple, analytic model for the enve- 
lope of the interference pattern observed in Fig. 01 Down- 
ward emission is generated in two pathways, namely from 
|2, fc — q) to |1, 0) and from |2, fe) to |1, q). We concentrate 
on the first pathway. To obtain a simple estimate for the 
envelope, we assume that each atom in state |2, k — q) 
will emit exactly one photon propagating downward if at 
the position of this atom there is a nonzero density of 
atoms in state |1, 0). This neglects the competition with 
rightward emission and the interference with the other 
pathway. 

Hence, we need to calculate the number Nj*™ 1 * of atoms 
in |2,fe — q) that have spatial overlap with any atomic 
density in 1 1 , 0) . For the following integration, we choose 
the coordinate origin such that the two relevant BECs lie 
on the y axis at positions ±y . iV™ 1 ^ is given by an inte- 
gral over the density n(x) of the BEC |2, fc — q) centered 
at y = —yo- We approximate all BECs as Thomas- Fermi 
parabolae with identical Thomas-Fermi radii. According 
to Eq. ([T]) we obtain in the new coordinate system 



n{x) 



Nk-q 

V 



x 

I?, 7. 



(y + yo) 
R l 



Rl) 



(25) 



where this is positive and zero otherwise. The integration 
boundaries have to be chosen such that only those atoms 
are counted that lie within the Thomas-Fermi parabola 
of the BEC 1 1,0) centered at y — y . For < y a < R y , 
the integration boundary is set by the BECs at y = yo 
and y = —yo for y < and y > 0, respectively. We 
substitute x = R x p cosy and z = R z psmip. We obtain 
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for < j/o < R y 



-Ry-ya r^i-(v+yo) 2 /R 2 y 
2nR x R z N k _ q J dy J dpP < x) 



3 ry/i-(y-yo) 2 /R 2 y 

dy / dp pn(x). (26) 

(Ry-yo) Jo 



This yields for all yo 



^k-q _ , ( yo_\ _ , I £fwm\ 



\ 2t J 



(27) 



where we used that the distance between the two relevant 
BECs amounts to 2yo = 2?ifc s iFWM/™ = -Rj^fwmAo 
with t Q from Eq. (|12l) . In addition, we introduced 

fc(a)= f (l- aJ )»(l + 3x + *») ) 0^<1 (28) 

and h(—x) = h{x). Within our approximation, the frac- 
tion h of the atoms in |2, k — q) will emit downward. A 
similar argument for the second pathway yields the same 
h. 

As Nk-q is given by an expression that is analog to 
Eq. (fl7|) . we obtain a simple estimate for the envelope 
of the downward emission by multiplying Eqs. (|TT)) and 
(1271) . to appears as a fit parameter in both equations. 
We decouple these two fit parameters by choosing the 
following fitting function for the data in Fig. [4] 



.4 



h(t) 



tr 



h [ ^-j (l + Fcoscj(i-t 5 )).(29) 



Here yo and t§ are offsets, A is an amplitude, and we 
expect f 4 = t . A fit of this model yields a curve that 
looks pretty similar to the curve shown in Fig. |U The 
best-fit values V = 0.38 and w/27r = 15.4 kHz agree well 
with the results from the simple fit shown in Fig. SJ With 
the best-fit values to = 0.97 ms and t$ — 0.42 ms, the 
envelope itself is well approximated by a Gaussian. We 
attribute the deviation from the prediction to — t^ — 
2.3 ms to the coarse approximations used to derive the 
envelope. 



III. EXPERIMENTAL DETAILS 

A. Rabi Frequencies 

The parameters in our experiment varied between 
different measurements. As an example, we list the 
parameters for Fig. |4l The applied laser powers are 
(P s ,i,F s , 2 ,P C 4,P c , 2 ,P c , 3 ) = (1,1,5,75,180) /xW. Here 
and in the following, the indices s and c refer to sig- 
nal and control light whereas the indices 1, 2, and 3 re- 
fer to the pulse number. The beam waists (1/e 2 radii 
of intensity) are (w Sl i, w Cl i, w c _2) = (0.17,0.32,1.8) mm. 
The use of identical beam paths yields w s ^\ = w St 2 



and w Ct \ — w Ci 3. From these parameters, one can es- 
timate intensities and electric field amplitudes according 
to / = 2P/ttw 2 and £ = ^2I/ce , where eo is the vacuum 
permittivity. We first consider the closed cycling transi- 
tion \F = 2,mj? = 2)H \F' = 3, m' F = 3) on the D 2 line 
at A = 780 nm. Its excited-state decay rate T = 1/(26 ns) 
yields a dipole matrix element d cyc = 2.5xl0 -29 Cm and 
a saturation intensity 7 sa t =1.6 mW/cm 2 . If the light 
drove this transition, then the Rabi frequency would be 
given by il^ = d^S/H = 

The detunings of the control light from the |2) — > |1') 
transition for the three pulses are (Ai,A2,As) = 2ir x 
(406,406,512) MHz. The effective two-photon Rabi fre- 
quency in our experiment is the superposition of the con- 
tributions from the two excited states |1') and |2') 



QcycQcyc 




V3 



12 A - A 



HFS 



(30) 



with the excited-state hyperfine splitting Ahfs/2tt = 812 
MHz on the D 1 line at A = 795 nm. The factors ±V3/12 
represent the products of the dipole matrix elements of 
the involved transitions in units of d cyc . 

With pulse durations of (ti,t 2 ) = (23,35) ps, we esti- 
mate pulse areas $ = Q c st of (^1,^2) = 7r x (0.49,0.51) 
which agree fairly well with the observed atomic popula- 
tion transfer. 

For retrieval, the parameter g from Ref. [S4j corre- 
sponds to n o g with the electric field amplitude £ s ,3 = 
yjTioj /2eoV q , the quantization volume V q , and uj — 2ttc/X. 
The connection between a homogeneous model with V q 
and a Thomas-Fermi parabola is made by approximating 
the quantization volume V q by the BEC volume V from 
Eq. ([21) • This is because the relevant physical quantity is 
the typical atomic density which equals N/V q and N/V, 
respectively. We obtain g = f2 e ff,3 = 2ir x 0.6 MHz. 



B. Raman AMW during the second Raman pulse 

During the second Raman pulse, Raman AMW occurs 
as a side effect with atoms transferred from |2, fe) to |1, 0) 
under emission of downward propagating light. The ef- 
fect on the atomic populations is negligible. But the 
emitted light reaches our EMCCD camera which needs 
~ 0.1 ms to erase this signal from the relevant pixels. 
This is why no data are shown for shorter times in Fig. 

SI 



C. No Raman AMW during the depletion pulse 

The depletion pulse propagates downward. After the 
corresponding absorption recoil, there is no emission di- 
rection available such that Raman AMW could occur. 
Hence, the depletion pulse causes only spontaneous light 
scattering. The depletion pulse depletes the F = 2 states 
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and operates on an open transition so that an apprecia- 
ble fraction of the atoms scatters only one photon. This 
does not suffice to spatially remove the atoms quickly 
from the BEC. Nevertheless, the spontaneous scattering 



aborts the FWM, because FWM occurs with BECs, but 
not with atoms incoherently distributed over a large num- 
ber of momentum states. 
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